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SUMMARY 


This papei' deals with the effect of the spectrum of the forcing fiuic- 
tion on the human pilot dynamics in manual control. A simple compensatory 
tracking experiment was conducted, where the controlled element was of a 
second-order dynamics and the forcing function was a random noise having a 
dominant frequency. The dominant frequency and the power of the forcing 
function were two variable parameters during our experiment. 

The results show that the human pilot describing functions are dependent 
not only on the dynamics of the controlled element, but also on the charac- 
teristics of the forcing function. This suggests that the human pilot be- 
havior should be expressed by the transfer function taking into consideration 
his ability to sense and predict the forcing function. 


SYMBOLS 


Aij(k) 

element of k-th autoregressive coefficient matrix 

B 

backward shift operator 

c(t), c(n) 

human pilot output 

dB 

decibel 

e(t), e(n) 

displayed ei*ror 

i(t), i(n) 

forcing function 

Kf 

static gain of forcing function filt<;r 
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M 


order of autoregressive model 
controlled element output 


m(t), m(n) 

s 

Yj(jw) 

Yf(Ju)) 

Yp(juj) 

A 

U)f 

Wn 


variable of Laplace transform 
controlled element 
forcing function filter 
human pilot describing function 
sampling interval 

damping of forcing function filter 
damping of controlled element 
power of forcing function 

undamped natural frequency of forcing function 
undamped natural frequency of controlled element 


INTRODUCTION 


It is well knoxm that when a hiiman pilot controls the system, his 
control behavior depends on the characteristics of the forcing function to 
the system as well as of the controlled element itself. A great number of 
papers have been published on this problem. 

Concerning the effect of the characteristics of the controlled el -nent 
on pilot behavior, Washizu and Miyajima (reference 1), and Goto and Lzu 
(reference 2 ) pointed out in a series of study on manual contrt>l of a second- 
order system that the human pilot takes notice of the periodicity in the 
response of the controlled element, if any, and makes use of it to improve 
his control performance. 

On the other hand, concerning the effect of the forcing function on 
pilot behavior, McRuer and Krendel (reference 3) pointed out that as the 
bandwidth of the forcing function increases, the effective time delay re- 
duces probably due to the muscular reaction characteristics of the human 
pilot. 

The purpose of the present paper is to investigate the effect of the 
forcing ftinction spectrum on the himian pilot dynamics in manual control. A 
simple compensatory tracking experiment was conducted, where the controlled 
element was of the second-order dynamics and the forcing function was a 
random noise having a dominant frequency. The dominant frequency and the 
power of the forcing function were two variable xjarameters during the experi- 


ment. Pilot describing 
coefficients identified 


functions were derived from the autoregressive model 
using the Akaike's Final Prediction Error method. 


experiment 


The system of our experiment was built up with an analogue computer, an 
oscilloscope and a control stick with a restoring spring. Its block diagi-am 
is as shown in figure 1. The error e(t) was displayed on the oscilloscope 
by a line segment moving vertically. The pilot was requested to minimize 
the error to the best of his ability. The controlled element had a second- 
order stable dynamics, and its transfer function was of the form; 


Yc(s) 



2 2 
s +2CnWns-H»)n 


( 1 ) 


The damping Cn and undamped natural frequency bin of controlled element 
were held fixed throughout our experiment such as, 

C = 0.1 

n * 

u) */20 = 4.1)7 (rad/sec) . 


The shaping filter of the forcing function also had a second-order stable 
dynamics as. 


Y^(s) 


^f^f 

2 2 


( 2 ) 


where the damping Cf was held fixed to 0.1 and the static gain Kf and the 
undamped natural frequency 03f. were two variable parameters. Thus, the white 
noise was transformed into a forcing function having a dominant frequency 
after passing the filter. The dominant frequency was varied by selecting the 

values of 0)^ as. 


0 )^ = 3.16, 2.24, 1.58 (rad/sec) . 

We chose four levels for the power of the forcing function by adjusting 
Kf of the equation. 
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( 3 ) 


0 ) K 
f f 


where is the power of the noise source. 

The expeiiment was of 12 cases, namely 3 kinds of frequencies and I 4 
power levels of the forcing function, and two runs of each case were perform- 
ed. After sufficient excercise, the analog data of the le»igth of 90 
seconds for each runs were recorded. The data, i(t), e(t), c(t) and m(t) 
in figure 1 were transformed into digital data by use of the NOVA mini- 
computer system. The FACOM 230-75 computer was employed for numerical calcu- 
lations of the following time series analysis. 


ANALYSIS 


By the use of the experimental data thus obtained, the human pilot 
describing functions were identified utilizing a time domain technique; tl.it 
is, an autoregressive model was fitted to the data by using the Akaike's 
MFPE (Multiple Final Prediction Error) method, (reference k) 

In the first place, the data were sampled from the analog data of the 
pilot output c(t) and the error e(t) with the sampling interval A, which was 
set as 0.1 sec. The sampled data are denoted by c(n) and e(n). Then, the 
autoregressive model of the form; 


✓ V 




✓ N 





c(n) 




c(n) 

+ 

°11 '^12 


C^(n) 

e(n) 


Api(B) A^jCB) 


e(n) 
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A^.j(B) = a^j(l)B + a^^(2)B^+ + a^j(M)B^ (5) 


Bx(n) = x(n-l) 


( 6 ) 


was fitted to the given data. B is the backward shift operator as shown in 
equation (6), and A£j(B)'s in equation (4) are the power series in B that 
are made up of the autoregressive model coefficients aij(k) with k going 
from 1 through M. The order of the model M is determined by the MFPE 
method. Ci(n)'s in equation (4) are mutually independent white noises. 

Once we have succeeded in fitting the model to the given data, namely, 
^12“ ^'’ 2.1 0) we can compute the pilot describing function using 
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/N A (jto) 

“ 1 - A^^(jw) 


( 7 ) 


where A 3 ^ 3 ^(jw) and Ai 2 (jw) are obtained from Ah(B) and Ai2^^^ equation (U) 
respectively, by replacing B with exp(-jwA). 

This method has recently been put into practical use, and our experience 
in using it has proved that it is quite efficient and powerful (reference 5)» 
Application of this method to our data was also suc cessful , as the estimated 
correlation coefficient of the noise sources, 0 ^ 2 / ^ ^11<^22> quite small. 


RESULTS 


Figures 2 and 3 are examples of the time histories of the records. 

Note that in figure 2, namely when the frequency of the forcing function (O^ 
was large, it is not evident that c(t) was affected by the forcing function 
periodicity. The pilot seemed to suppress only the controlled element 
periodicity. 

On the other hand, figure 3 shows the time history of the case when ih^ 
was relatively small. In this case, it is evident that c(t) was made up of 
two main sinusoidals ; one reflected the forcing function periodicity and the 
other reflected the pilot behavior which seemed to suppress the controlled 
element periodicity. This suggests that, when was relatively small, the 
human pilot behavior was affected by the forcing function. 

Above tendencies can be seen more obviously in the power spectri^ 
densities of the pilot output as shown in figure 4; namely in the vicinty of 
w = (Uf, the power spectra were pulled up as increased, and this phenome- 
non became more conspicuous when uip was relatively small. 

Typical pilot describing functions are shown in figures 5 and 6. From 
these figures, the following tendencies have been observed; 

1 ) If the power of the forcing function is increased, while keeping 

the undamped natural frequency unchanged, the gain of the pilot des- 
cribing function increases, but the phase lead becomes smaller in the 
frequency region below the undamped natural frequency of the controlled ^ 

element (%. g 

2) If the frequency of the forcing function Wf is decreased, while keeping | 

the power unchanged, the gain of the pilot describing tunction 

increases, but the phase lead becomes smaller in the low frequency range, 
especially in the neighbourhood of the undamped natural frequency of the 
forcing function. 
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Figure 7 shows the performance of the pilot control indicated by o2/o2. 
IS evident that the smaller t»ie undamped natural frequency was the ‘ 
better the performance became. This implies that when Wf was sLll/the 

becMc his 

^lese results lead to the following consideration concerning the for- 
cing function effects on human pilot control behavior. 

f-.-h effect of the forcing function bandwidth on the pilot describing 

fnr'r reference 3. It is pointed out in the report that 

forciJI''funct\'nr^ decreases as the bandwidth of the 

lorcing I unction increases* 

On the other hand, the present study put emphasis on tlie effect of the 

SS Sar?."*--"" »as beL 

to^wo^k o / increase in the power of the forcing function is likely 

the dominLr employ the control that takes into account 

periodicity in the forcing function. The attempt to suppress 

the f^fl^ency component may lead to the reduction of the power of 

suggested from the present study that if the 
response of the controlled element and the forcing function Lve perioLci- 

us^of^^h augment the system stability by making 

try to maL^tr" i response of the controlled element, and thfn, 

city of the Srcfnr?™M"-® possible by making use of the periodi- 

- ^ . . forcing function. Especially, if the power of the forcing 

function IS large and the two natural frequencies are separated, it would 

^ easy for the pilot to notice theee frequenclee and to'maje use of tSeL 
frequencies in the control, i-nese 

The present study has shown that the human pilot dascribing functions 
are dependent not only on the natural frequency of the controlled element 

'»■ forcing functloi TteT 
ults oeem to suggest that the human pilot control behavior couldn't be 

de?rr^^? function compensating the controlled element 

expressed by .the transfer function taking into 
consideration his ability to sense and predict the forcing function. 


CONCLUDING REMARKS 


such 1st of the forcing function on the human pilot 

1) If the power of the forcing function Oi^ increases, the gain of the 

i^an describing function |Yp| increases, but the phase lead of Y„ 
becomes smaller at o) < wn. P 

2) If the undamped natural frequency of the forcing function de- 

creases, the gam |Yp| increases but the phase lead of Y„ becomes 
smaller especially in the vicinity of ai = o)f. ^ 
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3) The human pilot seems to try to augment the system stability and 
make the performance better by use of and Wf, especially when 
is large, and and U)f are separated. ^ 
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diagram of* the Experiment 








Figure 2. An Example of the Time Histories of the Recorded Data 
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2,2^Krao/sec) 








